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It is now well known that electron (liole)-vibron coupling and hence Jahn- Teller (JT) effect is 
important understanding the properties of Ceo and related molecules. In this paper, we study 
H (g) {2h © g) coupling case to find the potential energy surfaces for the positively charged Ceo 
molecule due to distortion. The H 8 (2/i © g) Jahn- Teller system is of particular importance as 
this will be the JT effect displayed by Cgo molecules removed with an electron. Cgjj is obtained by 
removing one electron from fivefold degenerate Hu highest occupied molecular orbital (HOMO) and 
a hole in HOMO interacts with the vibrational modes of Ceo and symmetry is broken. We apply the 
method of symmetry breaking mechanism to obtain expressions for the potential energy surface. 

PACS numbers: 

INTRODUCTION 

Except at absolute zero in temperature, all atoms and molecules have energy and vibrate. The electrons in Ceo 
molecule will be sensitive to these vibrations. It is believed that this vibronic coupling plays an important role in 
determining the behaviour of Cgo molecules and related fuUerene compounds. 

The symmetry group of Cqq molecule is icosahedral (Ih) ■ This symmetry group implies large representations, 
thus large degeneracies of the interacting electronic and vibrational states of isolated ion. Icosahedral symmetry is 
extremely rare in nature, vibronic coupling effects have not been investigated in this symmetry until very recently. 
Many interesting effects due to vibronic coupling are possible from a theoretical point of view due to quantum- 
mechanical four and five-fold degenerate states. The ground state of pure Ceo molecule is singlet Ag state, so it is not 
very sensitive to JT effect. However, the ground state of C^q and C^^ states do potentially strong effects. Ionic form of 
Cgo molecule is subject to the Dynamical Jahn- Teller effect (D JT) . The Dynamical Jahn- Teller effect is closely related 
to, but should be distinguished from the static Jahn- Teller effect, where there is a permanent symmetry breaking 
molecular distortion or instability to one of the potential-well minima, removing all or part of orbital degeneracy. 
JT effect is also important from experimental point of view. Most workers agree that electron-vibron interactions 
are very important for many properties of fuUerenes and may explain superconductivity. It is known that origin of 
superconducting pairing in Cgo compounds is due to interaction of electrons with vibronic modes of Ceo molecule. 

In the general formalism of the JT effect, a degenerate electronic state corresponding to a representation of D of 
symmetry group G of the molecule can interact with the vibrational modes corresponding to representations contained 
in the symmetric part of the direct product D(i)D (excluding the identical representation which is trivial). As it is well 
known, the molecular symmetry reduced by the JT distortion with splitting of the electronic-state degeneracy. The 
distorting forces acting along a certain non-totally vibrational modes carry the nuclei over into distorted configurations. 
Distorted configuration of a molecule can be characterised by subgroup symmetries of parent molecular group. 

In general, a physical system may pass from a symmetric to less symmetric state, during its evolution. This 
symmetry transition is known as spontaneously symmetry breaking [6,7]. The process of symmetry breaking applies 
to all domains described by expectation values of operators. This common mathematical technique can be applied 
successfully to different domains of physics because they all share in the occurrence of broken symmetry. 

In order to study the JT structure of A^Cqo and/or C^" various theoretical approaches have been used. Ceuleman 
has proposed on analytical treatment of JT distortion in the general case of a fivefold degenerate state of an Ih molecule. 
The extreme points were identified by using the isostationary function and epikernel methods [1,2]. First order JT 
interaction and its continuous group invariance were discussed and the energies were found for JT coupling with 
and without splitting [3] . The ground state Berry phase of some JT systems were calculated in the five dimensional 
rotation group SO (5) [4]. Classical and semi-classical models were introduced about JT manifold for electron vibron 
interactions, by Auerback, Manini and Tosatti [5]. The H ® {h(B g) problem is investigated analytically using a unitary 
transformation by CP Moate, J.L Dunn et al. [12]. 

The present paper provides a group theoretical treatment of JT distortion in general case of a fivefold degenerate 
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state of Ifi molecule. The JT surfaee energies have been obtained by breaking symmetries of Ih group into its maximal 
little group. This is never done before. 

SYMMETRY ASPECTS AND COUPLING OF STATES 

In this section, we start by describing the Hamiltonian that generates DD surface where D denotes the irreducible 
representation. The standard Hamiltonian may be written in the form 

H = Ho + Hjt; Ho = lh^{Pi + Q^) (1) 

i 

where Hq describes free (uncoupled) electrons/holes and vibrations and Qi is the distortion coordinate and P, is 
the conjugate momentum. In general, Hjt introduces a rotationally invariant linear coupling between the electronic 
state and the vibrational mode. It is known that the Hamiltonian for a linear JT coupling oi H ^ (2h ® g) system is 
invariant under the rotational operations of 1^ group. If either the restriction of linear coupling or harmonic forces is 
relaxed, the symmetry of Hamiltonian is reduces to point group G. 

As we have mentioned above, totally symmetric part of direct product of an irreducible representation of a finite 
group, which describes the properties of JT surfaces, is written in the form of 

[£>^ O £>^] = ® £>^^ ® • • • ® (2) 

where is the angular momentum quantum number. Decomposition of [D^ (g) D^] implies that the JT Hamiltonian can 
be written in the following way 

Hjt = H^' +H^' + ■■■ + H^- (3) 

where is the JT Hamiltonian and it is invariant under the symmetry operations of corresponding finite group, for 
the 2^+1 dimensional representation. Since the Ih group is subgroup of 0(3), decomposition of coupling of two states 
should be written in terms of 

Therefore, symmetric part of the five dimensional direct product for Ih group can be written in the notation of 
icosahedral molecule: 

[H^ ® = [Hi] = [Hi] =Ag®Hg® {Gg 9 Hg) (5) 

where Ag,Gg and Hg are one, four and five dimensional irreducible representations of Ih, respectively. The decom- 
position of direct product of ig) Hy] coupling in 0(3) is found using the Table 1 and is in the form 

HjT = H'' + H^ + H^ (6) 

where the superscripts 2 and 4 are £ values. The symmetric part contains the totally symmetric representation 

Hq = Ag that is trivial polaronic problem [9,10] and it can be solved, exactly. The Hq vibration shift energies 
but it does not cause splitting. Remaining vibrations subtend the configuration space, which contains all distorted 
configurations, may be reached by JT active coordinates. As shown from Table 1, the Hamiltonian H2 corresponds 
five dimensional representations Hg and H4 corresponds to the direct sum of GgHg vibrational levels. 

DETERMINATION OF STATIONARY POINTS ON JT SURFACE USING LITTLE GROUPS OF 

ICOSAHEDRAL 

The second important aspect of invariant polynomial functions, in addition to, linear Jahn- Teller matrices, concerns 
the extremum points on JT surface. Since our interest are the minimal points and the saddle points, the present 
section is devoted to the determining of extremum points on JT surface by breaking symmetries of icosahedral group 
into its maximal little groups. The problem has been studied by Ceulemans by using the isostationary function [1]. 
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Ih 




Ih 





Ao 




7 Tiu + T2u + Gu + Hu 


1 


Tiu 




8 T2g+Gg + Hg+Hg 


2 


H, 




9 Tiu + T2u + Gu + Gu + Hu 


3 


Gu 


+ T2u 


10 Ag + Tig + T2g +Gg+ 2Hg 


4 


Go 


+ H, 


11 2T,„ + T2u + Gu + 2Hu 


5 




+ T2u + Hu 


12 Ag + Tig + T2g + 2Gg + 2H g 


6 A, 


+ Tig + Gg +Hg 


13 Tiu + 2T2u + 2Gu + 2Hu 



TABLE I: The representations of the 0(3) group with i <14 are spUt by the icosahedral point group. 



D2 D3 D5 


T 


H 

G 

H®G 


2A + B1+B2 + B3 Ai +2E A1 + E1+ E2 

A + B1+B2 + B3 A1+A2 + E Ei + E2 

SA + 2Bi + 2B2 + 2B3 2Ai + A2 + 3EA1 + 2Ei + 2E2 


E + T 

2A + E 

2A + 2E + T 



TABLE IL Maximal Little groups of H and G representations. Decomposition impUes that T is not maximal little group of H 
state, and Ds is not maximal little group for G state. 



The little group is the new symmetry group of distorted molecule as a result of coupling. A real representation of 
any subgroup S < G; the degree of subduction can be computed by the relation 

pes 

where Xz{p) is the character of the representation p. In a 2;th representation of finite group, if all subgroups S > S and 
C(S) < C{S) then S is little group of G. For finite groups converse of this definition is also true. In order to simplify 
the notation, we shall switch from 0(3) to the SO{3) notation, omitting therefore the g/u symbol for inversion. Then, 
the electronic-vibrational coupling is symbolized as {2h(Bg) where capital letter shows electronic and small letters 
show the vibrational levels. Without loss of generality the problem for Ih group can be treated equally well in the 
subgroup of proper rotations of group /. The maximal little subgroups of /group have computed by using Eqn.(7) 
and are given in Table 2. As seen from Table 2, the maximal little groups predict the existence of dihedral groups 
1)5,1)3, and D2 minima on the icosahedral molecular cage. In the coordinate space the distortion that is the result 
of JT instability should conserve all maximal little group symmetry. 

In our perspective the structure of Jahn- Teller surfaces have been identified by the symmetry breaking of the 
continuous symmetry to the true finite point group of the representation space and its maximal little groups. This 
may be represented as follows 

I^S' (8) 

where S' are little groups of / group given in Table 2. 



GROUP INVARIANCE OF JAHN-TELLER SYSTEMS 

As we stated in the previous section, the Hamiltonian of (2/i ® g) coupling have three parts H'^ and H'^, 
which must be separately invariant under / symmetry. In this section, we want to construct a polynomial function in 
electronic and nuclear configuration space, to examine symmetry properties of potential energy surface. A polynomial 
function which have been produced for this purpose is in the form 

2i+l 2m+l 

UUX,Q)=J2 ^iokXiXjQk. (9) 

— l k—1 

In this expression Xi and Xj correspond to electronic coordinates Qk corresponds to nuclear coordinates. The force 
elements or coupling coefiicients should be chosen appropriately for each H ® h and H ® [h® g) coupling. In the 
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Eqn.(9) the indices 2£+l and 2m + 1 stand for dimensions electronic and active coordinates, respectively. For H ®h 
coupling f and m take value 2 and for H ® (h® g) case£ and m take values 2 and 4 respectively. 

We will focus our treatment for the computation of invariant polynomials for iJ (g)/i and H^{h(Bg) problems. The 
invariant polynomial function for £ = 2, in real basis has been computed using the matrix representations of H state 
which are given in Appendix. The 5x5 generators transform electronic coordinates Xi(i = 1, • • • ,5) and nuclear 
coordinates Qj. It can be written in the form 



(10) 



i,j=l k—1 



5 5 

In this equation, Xj = ^ ^In-^n and Q'k = J2 ^InQn- is the matrix elements of 5 x 5 generators of /. The 

71=1 n=l 

Eqn(lO) is solved for Fijk and two linearly independent polynomial function have been obtained. One of our main 
goal that first order JT interaction matrices can be derived by working out invariant polynomial function. The double 
differentiation of U22iX,Q) with respect to electronic coordinates Xi and Xj produce linear JT interaction matrix. 
In general, we can write 



d^U2m 

dXidXj 



The JT interaction matrix for (g) /i coupling is 



So = i^i 



/ 2Q5 


-7303 

\ 2Qi 
/ Vl5Qi- 





-V5Q3 
2Q2 

#04 4 



-V3Q3 

-VSQi - C 
-V3O2 



#05 



I5O2 + #03 

— O2 



V304 

-\/303 

-V3Q2 

V3Qi - O5 
-O4 

T5O2 + #03 



= 1,2- ■ -5). 



2O1 \ 
2O2 

-O3 + 
-O4 
-2O5/ 



(11) 



(12) 



V 



-#0i 
#0i- 



-305 
304 



/2 -I- -^W3 
-VT5Q1 + #04 + #05 
#0i - 305 

— W2 



#02 



#02 - 303 



#0i - 305 
-2^/304 - \/505 

-2V^Q3 
-3O2 - V503 



■#0i 



305 



— O2 
-2\/303 

2^304 - 7505 
-30i - 7504 



The force elements Fi and F2 are given by 

1 



Fi 



135^3 



{-ISVEFhi + 15Fh2); F2 = 



1 



27V15 



{VEFhi + 3Fh2) 



#0i - 304 

#02 - 303 
-3O2 - 7503 
-30i - 7504 

7505 



(13) 



In this equation F^i and Fh2 are coupling parameters of H mode. The invariant polynomial function is derived 
for H {h(B g) as in the same way H ®h. In this case Eqn.(lO) takes form 
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^FijkXiXjQk = ^ FijkXlX'jQ'^. 



i,j=l k=l 



(14) 



i,j=l fe=l 



Nine dimensional H iSi {h(B g) state consists of one H mode with components {Oi' O2. Q:h Qi, Os} and G mode with 

components {Oe; Oti Os) Og}- The icosahedral generators transform the electronic and nuclear coordinates. Trans- 
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formation Xi and Xj are same as given in Eqn.(lO). The nuclear coordinates Ofc is transformed as Ofe = ^L^"' 

n=l 

where is the direct sum of 5 and 4 dimensional irreducible matrix generators T and A given in Appendix, respec- 
tively. Solution of the Eqn.(14) for coefficients Fijk gives that three linearly independent function. In Ceulemans's 
paper [1], relations between polynomial coefficients Fijk in Eqn(lO) and Eqn(14) arc expressed in terms of Clebsch 
Gordon series for the icosahedral point group [13]. The linear JT interaction matrices for this coupling are derived 
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from the relation (11). The two of them are same with matrices oi H h coupling and third one is 



F3 



( -Qe-2V5Qs 

2V5Q7 



V 



IQ7 - "^Qo 



2V5Q7 

-Qe + 2V5Qs 



jQr — ^^Qg 

4Qe + V5Q8 
V5Q7 



IQ7 + ^Qg 

V5Q7 
4Q6 - V5Q8 
-¥Q8 



-6Q6 / 



(15) 



The sum of the B2 given in Eqn.(12) and JT interaction matrix given in Eqn.(15) corresponds to the first order 
JT interaction matrix (B4) for H ^ {h (B g) coupling. It is obvious that the interaction matrix is also obtained by 
considering only H ® g coupling. Force clement ^3 is related by coupling parameter of G mode and is given by 
-P3 = -Pg/9- The first order linear JT interaction matrices for H state have been found in [14] and are in agreement 
with our results. We guess that the higher order JT interaction matrices may be obtained by constructing higher 
order icosahedral invariant polynomials. 

In section 4.1, symmetry of icosahedral is broken into its little groups iov H ® h and H ® [h® g). Combination of 
eigenvalues of the with harmonic potential energy, in terms of little groups yield JT surface energy. 



TRANSITIONS ASSOCIATED WITH H^h AND H^{h®g) COUPLINGS 

The five dimensional irreducible representation of I group has three maximal little group named D2, D3 and D5. 
In order to break symmetry of a parent group into its little groups, one should assign an appropriate Qi which can 
be computed by constructing set of equations such that 

2e+i 

where V^j is the matrix elements of generator of the corresponding little group. The method of symmetry breaking 
predicts the existence of saddle points, trigonal and pentagonal turning points on JT surfaces associated with 
and groups. From the solution of Eqn.(16) it is found that two H type coordinates QhiiQh2 and one G type 
coordinate Qg- The H and G type coordinates are constrained to Q4, = —tQhi, Q5 = —-^Qh2, Qe = Qg, where 
r = 5(1 + \/5). Under the given conditions; the energy eigenvalues of each little group are computed. 

4.1.1 D2 TRANSITION 

Decomposition of five dimensional representations in D2 group is 2 A + Bi + B2 + B3. The symmetry of the group / 
is broken into 1)2, assigning as Qi —> a/5(34 — \/3(35, Q2 ^ 0, Q3 ^ 0, Qa^Qa, Q5 ^ (55, using the Eqn.(16). 
After substituting values of Qi into B2, the eigenvalues of B2 arc carried out. Combinations of eigenvalues of B2 with 
harmonic restoring potentials for the distortional coordinates Qhi, and Qh2 gives JT surface energy values. 

E{A) = ±^ [(3i^li + Fl2){2T^Qli - 5tQhiQh2 + 5QI2)] ^ + ^KniQl, + QI2) 



E{Bi) = -^{Fhi + Fh2){^tQhi - 5Qh2) + (Fhi - ^Fh2)V5Qh2 + \kh{Qhi + 

E{B2) = -^{Fh2 - F„i){4tQhi - 5Qh2) + {Fhi + 3Fh2)V5Qh2 + ^^^^(Qhi + Qh2) 
4V5 2 

E{B3) = -^{2TFHlQHl-l{FHl-V5FH2)QH2 + ^KHiQli+Ql2)- (17) 

In this equation, Kh is the harmonic force constant. The energy values of E{A) predicts the existence of saddle 
points. 

The direct sum of g and h states consists of both G-type and i?-type nuclear coordinates. For this reason, 

computations are more complicated than the h state. In this case, there arc three distortional coordinates that are 
invariant under D2. In Eqn.(16) nine dimensional matrix generators which have been obtained by taking direct sum 



6 



of five and four dimensional representations are used to assign Qi values. It is found that Qi \/5Qi — V^Q^, Q2 
0, Qs ^0, Qi, Q5 Q5, Qe Qe, Q7 0, Qs ^ V^Qe, Q9 0. In this basis, combination of eigenvalues 

of B4 with harmonic potential gives us 

E{A) = ±^ [{3F^, + F^^){2t'qI, - 5tQhiQh2 + 5Ql^)] ' + ^FgQg + ^KcQl + \kh{.Q\^ + QI2) 

E{Bi) = -^{Fhi + Fh2){^tQhi - 5Qh2) + {Fhi - 3Fh2)\/5Qh2 - \fgQg + \kgQI + \kh{QIi + QI2) 

E{B2) = -^i.FH2 - Fhi){AtQhi - 5Qh2) + {Fhi + 2,Fh2)VIQh2 - \fgQg + \kgQg + \kh{Q\i + QI2) 

E{Bs) = -^{Fhi22Qhi - ^(^hi - VEFh2)Qh2 - \fgQg + \kgQI + \kh{Q%x + QL), (18) 

where Kg in the Eqn.(18) is the harmonic force constant of G mode. It is obvious that the magnitude of splitting of 
energy values is increasing in hg state compared to h state. 



4.1.2 D3 TRANSITION 



In this case decomposition of H representation of irreducible representations of D3 is Ai + 2E in D3. Following 
the same procedure given in 4.1.1, symmetry of / group can be broken into Ds, in the directions, Qi 0, Q2 ^ 
0) Qa ~* 0, Q4 ^ 0, Q5 — > Q5. The computed energy values are given as 

E{A,) = -^FhiQh2 + IkhQ%2 

E^{E) = \Fh,Qh2 T ^(^)'+(^)'Qh2 + \KHQ]t2- (19) 

Transformation of nuclear coordinates for H ® [h® g) state gives that Q\ ^ 0, Q2 ^ 0, Q3 ^0, Q4 ^ 0, Qs 
Q5, Qe — * Qe, Qi — * 0, Qs — * 0, Qg — > 0. The energy values consists of the two distortional coordinates, and 
therefore 

E{A^) = -1{FgQg + FhiQh2) + \kgQI + \khQ%2 

+ ]^KgQI + \khQ12- (20) 
These energy values are in agreement with the energy values given in [1]. 



E^{E) 



-^{FgQg + FhiQh2)t 

D 



IfhiQh2 - ^FgQg) + \Fl2QH2 



4.1.3. D5 TRANSITION 

In the quintet irreducible representation state reduces to Ai + Ei + E2. Applying the same procedure as in 
the previous section, Qi values are found as Qi — + — Q2 — * 0, Q3 — + 0, Qa 0, Q5 — » ^Qi- The 
corresponding energies are obtained as 

E{A^) = -1=Fh2Qhi + \khQI^ 

E{Er) = ^{-Fhi + ^Fh2)Qhi + \khQIi 

E{E2) = ^{Fhi + ^Fh2)Qhi + \khQIi- (21) 

The calculations have been carried out for H ^{h® g) coupling and same results have been obtained as expected. 

Tetrahedral group (T) is also maximal little group of G representation of I group. In H ® {h® g) coupling some 
energy values have been expected. Since (T is not little group of H representation. Thus va. H ® h coupling energy 
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eigenvalues of B2 is zero.) Decomposition of H state in T group is E + T, and symmetry is broken according to the 
directions Qg ^\f^Qei Qi (i = 1, . . .8) . The corresponding eigenvalues of B4are FqQg and — |JgOg, 

for E and T respectively. 

In general, minimal energy values for each little group can be obtained from the given energy expressions. 



CONCLUSION 



In summary, we have shown how the symmetry breaking method is applied for the determination of the potential 
energies of the {2h(Bg) surface. In Ceulemans' paper [1], these energies were found by the method of isostationary 
function and potential energies of and D5 groTips were investigated. In our work, all maximal little groups of 
icosahedral group arc studied for the H ^ (h Q) g) and H ® h state. Splitting of energy levels of icosahedral system 
due to distortion is analysed and amazingly interesting that, there is a proper contribution on the connection between 
our method and method of isostationary function. This method can also be used for other distorted systems. 
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APPENDIX 



Generators of /, D2,D^, D5 and T groups for Gand H representations 

Four and five dimensional irreducible matrix generators of / and its little groups can be generated from the matrices: 



r2 

p4 



I, 
D2 



^ T 

/ -2 




\ 3 



1^2,1^4 p2 p4 ^ jj^ 



V5 0\ 
3 
2 
0/ 



(2 






\ 
-1 \/3 

-V3 -1 
2 / 



( -4 a/is -\/5 \ 

-715 -3 \/3 -3 
-V5 -V3 1 3^3 
-3 -3^5 / 



/ 1 \ 
0-100 
10 

V -1 / 



6\/3 



/ 73 -3^5 715 -3^5 \ 

-3a/5 -6 -3a/3 

VIS -6 -4\/3 3 

V -3\/5 -3\/3 3 -3\/3 / 
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-2\/3\ 







3^5 


-6 










-6 


-3V5 







2V5 





-1 


2VT5 


V 


-2\/3 





2\/l5 


3 / 





-1 


-a/3 





o\ 




V3 


-1 
















-1 


V3 










-%/3 


-1 


V 











02/ 



36 



/ 7-9^5 
7\/3 + 3\/T5 

-6\/3 + 2\/T5 
18 + 2V5 
V 4V3 



-7V3-3VT5 6V3-2V15 I8 + 2V5 4V3 \ 
-21 + 3\/5 -6-6\/5 -6\/3 + 2\/l5 12 

-6-6\/5 3-3v^ -\/3-3yT5 -12\/5 
6^3-2^ V3 + 3VT5 -1 + 9V5 -4^/15 
-12 12^5 -4^15 12 / 
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/ 1 \ 

0-1 
-10 
1 

V 1 / 



12^3 



7 + ^/5 
IO + 2V5 



7 + 75 10 + 2\/5 -2^1 

-^/3 + ^/T5 -2V3-2^/T5 -14 + 2x/5 
-2V3-2VT5 7V3-VT5 -I-V5 
-1 - Vs 



2\/3 -14 + 2\/5 
V 2-6^5 -273-2^15 -2V3 + 2^/T5 
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-6-2^5 



2\/3 2 - 6\/5 \ 
2V3- 2VT5 
2V3 + 2VT5 
-6-2V5 

-4V3 y 
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